Ising model on nonorientable surfaces: 
Exact solution for the Mobius strip and the Klein bottle 
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Closed-form expressions are obtained for the partition function of the Ising model on an M x TV 
simple-quartic lattice embedded on a Mobius strip and a Klein bottle for finite M and N. The 
finite-size effects at criticality are analyzed and compared with those under cylindrical and toroidal 
boundary conditions. Our analysis confirms that the central charge is c = 1/2. 
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I. INTRODUCTION 

There has been considerable recent interest in 
studying lattice models on nonorientable surfaces, both 
as new challenging unsolved lattice-statistical problems 
and as a realization and testing of predictions of the con- 
formal field theory Q. In a recent paper |jj we have 
presented the solution of dimers on the Mobius strip and 
Klein bottle and studied its finite-size corrections. In this 
paper we consider the Ising model. 

The Ising model in two dimensions was first solved by 
Onsager in 1944 || who obtained the close-form expres- 
sion of the partition function for a simple-quartic A4xJ\f 
lattice wrapped on a cylinder. The exact solution for an 
M.x J\[ lattice on a torus, namely, with periodic bound- 
ary condition in both directions, was obtained by Kauf- 
man 4 years later || . Onsager and Kaufman used spinor 
analysis to derive the solutions, and the solution under 
the cylindrical boundary condition was rederived later by 
McCoy and Wu § using the method of dimers. As far 
as we know, these are the only known solutions of the 
two-dimensional Ising model on finite lattices. Here, us- 
ing the method of dimers, we derive exact expressions for 
the partition function of the Ising model on finite Mobius 
strips and Klein bottles. As we shall see, as a conse- 
quence of the Mobius topology, the solution assumes a 
form which depends on whether the width of the lattice 
is even or odd. However, all solutions yield the same bulk 
free energy. We also present results of finite-size analyses 
for corrections to the bulk solution, and compare with 
those deduced under other boundary conditions. Our 
explicit calculations confirm that the central charge is 
c= 1/2. 



II. THE 2M x N MOBIUS STRIP 

To begin with, we consider a 2M x N simple-quartic 
Ising lattice C embedded on a Mobius strip, where M, N 
are integers and 2M is the width of the strip. The ex- 



ample of a lattice C for 2M = 4, N = 5 is shown in Fig. 
1. 

While we shall consider the case of a uniform reduced 
interaction K, to facilitate considerations it is convenient 
to let the N vertical edges located in the middle of the 
strip to take on a different interaction K\ as shown. By 
setting K\ — the Mobius strip reduces to an M x 2N 
strip with a "cylindrical" boundary condition, namely, 
periodic in one direction and free in the other, for which 
the partition function has been evaluated by McCoy and 
Wu M . By setting K\ — oo the two center rows of spins 
coalesce into a single row with an (additive) interaction 
which in this case is 2K . These are two key elements of 
our consideration. 
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FIG. 1. A 4 x 5 Mobius strip C. Vertices labeled A, B, C, D 
are repeated sites. 

Following standard procedures we write the parti- 
tion function of the Ising model on C as 

Z^ N {K,Kt) = 2 2MN {co^K)^ 2MN - N \co S hK l ) N 
xG(z,zi), (1) 



where z = tannic", z\ = tanhiTi, and 



G(z,Zi) 



E 

closed polygons 



(2) 



is the generating function of all closed polygonal graphs 
on C with edge weights z and z\. Here, n is the number 



1 



of polygon edges with weight z and n\ the number of 
edges with weight z\. 

The generating function G(z, zi) is a multinomial in z 
and z\ and, due to the Mobius topology, the integer m 
can take on any value in {0, N}. Thus, we have 



JV 



(3) 



ni=0 



where T m (z) are polynomials in z with strictly positive 
coefficients. 

To evaluate G{z 1 z{), we again follow the usual pro- 
cedure of mapping polygonal configurations on C onto 
dimer configurations on a dimer lattice Co of 8MN sites, 
constructed by expanding each site of C into a "city" of 
4 sites [fi|f§. The resulting Cd for the 4 x 5 C is shown 
in Fig. 2. 




FIG. 2. The dimer lattice Cd corresponding to the 4x5 
Mobius strip. 

As the deletion of all z\ edges reduces the lattice to one 
with a cylindrical boundary condition solved in 0, we 
orient all edges of weights z and 1 as in M . In addition, 
all z\ edges arc oriented in the direction as shown in Fig. 
2. Then, we have the following: 

Theorem: 

Let A be the 8MN x 8MN antisymmetric determinant 
defined by the lattice edge orientation shown in Fig. 2, 
and let 



PL4(z,zi) = y/&etA(z,z{) 
denote the Pfaffian of A. Then 



N 



(4) 



(5) 



ni=0 



where e 4m = e 4m+1 = 1, e 4 m+2 = &im+3 = -1 for any 
integer m > 0. 



Remark: Define 

[N/4] 

X P =J2 T im+P (z)zt m+P , p = 0, 1, 2, 3, (6) 

m=0 

where [N/4] is the integral part of N/4 so that G(z, Zi) = 
X Q + X 1 +X 2 + X 3 . It then follows from (g) that we have 



PL4(z, ±i Zl ) =X + X 2 ± i(X 1 + X 3 ) 
As a consequence, we obtain 

i r. 



(J) 



G{z,zi) 



(l-i)PfA(z,izi) 
+(1 + i)PfA(z, -izi) 



(8) 



where, as evaluated in the next section, the Pfaffian is 
given by 



PL4(z,zx) = [z(l-z 2 )} 

N 

><n 



2\i MN 



n=l 



sinh(M + l)t(<f> n ) - c(z,Zx)smhMt(<f) n ) 
sinht {4> n ) 



, (9) 



with 



z(l + z 2 + 2zcos 4> n ) + 2(-l)"zi sine 



c{z,zi) - ! _ . , 

cosh t (4>) = cosh 2K coth 2K ~ cos q 
<j) n = (2n - l)n/2N. 



(10) 



Here we have used the fact that IIn=i = rin=ivy-i m ^ ne 
product in (||). Substituting these results into (Q), we are 
led to the following explicit expression for the partition 
function, 

1 / \ MN 
Z™° h N (K,K) = -(2 S mh2K) 

x [(1 - i)F+ + (1 + 



(11) 



where 



N 

F ± = il 

n=l 



/ e*^») -c(z,±iz) N 



_ e -Aft(0„) , 



2sinht((/)„) / 

e -t{<t> n ) - c (z,±iz) 
2smht(4> n ) 



(12) 



This completes the evaluation of the Ising partition func- 
tion for the 2M x N Mobius strip. Note that we have 
coshi(0„) > 1 so we can always take t(cf> n ) > 0. The 
leading contribution in ([l2|) for large M is therefore 

]-J e Mi(<A n )_ 

For the 2x5 Mobius strip, for example, we find 

PfA(z, zi) = l + z 10 + lOziz 5 - 5zfz 2 (l + z 2 + z 4 + z 6 ) 
-20zfz r ° + 5zfz 4 (l + z 2 ) + 2zfz 5 , 
G(z, zi) = 1 + z 10 + 10z lZ 5 + hz\z 2 {\ + z 2 + z 4 + z 6 ) 
+20 Zl 3 z 5 + 5zfz 4 (l + z 2 ) + 2zlz 5 , (13) 



2 



which can be verified by explicit enumerations. 
We next prove the theorem. 

Considered as a multinomial in z and z± , there exists a 
one-one correspondence between terms in the dimer gen- 
crating function G{z,z\) and (combinations of) terms in 
the Pfaffian . However, while all terms in G(z, z\) are 
positive, terms in the Pfaffian do not necessarily possess 
the same sign. The crux of the matter is to find an appro- 
priate linear combination of Pfaffians to yield the desired 
G{z, z\). For this purpose it is convenient to compare an 
arbitrary term G\ in the Pfaffian with a standard one Co- 
We choose Co to be one in which no z and z\ dimers are 
present. 

The superposition of two dimer configurations repre- 
sented by Co and C\ produces superposition polygons. 
Kasteleyn || has shown that the two terms will have 
the same sign if all superposition polygons are oriented 
"clockwise-odd" , namely, there is an odd number of edges 
oriented in the clockwise direction. 

Now since all z and 1 edges of Cd are oriented as in 
0, terms in the Pfaffian with no z\ edges (n\ — 0) will 
have the same sign as Cq. To determine the sign of a 
term when z\ edges are present, we associate a + sign to 
each clockwise-odd superposition polygon and a — sign 
to each clockwise-even superposition polygon. Then the 
sign of Ci relative to Co is the product of the signs of all 
superposition polygons. The following elementary facts 
can be readily verified: 

i) Deformations of the borders of a superposition poly- 
gon always change mi, the number of its Z\ edges, by 
multiples of 2. 

ii) The sign of a superposition polygon is reversed un- 
der border deformations which change mi by 2. 

iii) Superposition polygons having or 1 z\ edges have 
a sign +. 

iv) There can be at most one superposition polygon 
having an odd number of z\ edges (a property unique to 
nonorientable surfaces). 

Let mi = 4m + p, where m is an integer and p = 
0,1,2,3. Because of iv), we need only to consider the 
presence of at most one polygon having p = 1 or 3. It now 
follows from i) and iii) that = t± m +\ = +, and from 
i), ii) and iii) that en m +2 — £4m+3 = — • This establishes 
the theorem. 



III. EVALUATION OF THE PFAFFIAN 

In this section we derive the expression (^). 

From the edge orientation of Cd of Fig. 2, one finds 
that the 8MN x 8MN antisymmetric matrix A assumes 
the form 

A(z, zi) = A Q (z) <g> I 2N + A + (z) ® J 2N + A_(z) <8> j£ N 
+A 1 (z 1 )(g>H 2N , (14) 



where Aq, A+, A_, A\ are 4M x 4M matrices, I 2 n is the 
2N x 2N identity matrix, and J 2 n, H2N are the 2N x 2N 
matrices 
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(15) 



In addition, one has 

A (z) = a , O Im + a .i(z) ® F M + 00,-1(2) ® 
A±(z) = a±i >0 (z) <x> Im 
Ai(zi) = a{zi) (g) G M , 
where Fm , Cm are M x M matrices 



Fm — 



(16) 
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Gm — 



is the transpose of Fm, and 



«o,o — 



(17) 



a(zi) 



aoA z ) 
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a-i,o{z) 
00,-1(2) 



(18) 



ai,o{z) 
a o.i(z)- 

We use the fact that the determinant in (|j) is equal to 
the product of the eigenvalues of the matrix A. To eval- 
uate the latter, we note that J 2 n, J 2 n an( i H2N mutu- 
ally commute so that they can be diagonalized simultane- 
ously. This leads to the respective eigenvalues e l< ™ , e - *^™ 
and i(— l) n+1 and the expression 
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2.V 



det A(z,z{) = Y[detA M (z,zi',(j>n), 



(19) 



n=l 



where 

A M (z, z i; A Q (z) + A + (z)e^ + A_(z)e~^ 

+z(-l) n+1 A 1 (z 1 ) (20) 

is a 4Af x AM matrix. Writing out explicitly, we have 

A M (z,Zi;(t> n ) = 

( B(z) a ,i(z) \ 
a -i(z) B(z) 00,1(2) 





V 



a -i(z) B(z) 00,1(2) 
00,-1(2) C(z,Zi)J 

(21) 



where C(z, Zi) = B(z) + i{-l) n+1 a{z x ) , and 

-1 

_n _i_ „t,—i<t>n\ n 

B(z) = 



/ 1 + ze 1 ^ 

-(1 + ze"^™) 1 

1 -10 1 

V i 1-10 



(22) 



The evaluation of det Am (z, z%; tf> n ) can be carried 
out by using a recursion procedures introduced in |l(]] 
for a self-dual Ising model. Specifically, let Bm — 
detAM(z, Z\ \ 4>„) and Dm the determinant of the matrix 
Am(z, z\\ 4> n ) with the fourth row and fourth column re- 
moved. Then by expanding the determinants one finds 



the recursion relation (which is the same as that in |10|) 

M > 2, (23) 



Bm 



on 0,12 \ ( Bm-i 
a-21 a 2 2 J \ D M -i 



with 



an = 1 + z 2 — 2zcoscj) n , 
a\2 = — 2iz 3 sin^ n , 
ai\ — 2iz sin (f> n , 

a 2 2 = z 2 (l + z 2 + 2zcos0„), (24) 
and the initial condition (which is different from |hJ) 

Bi=B 1 (z,z 1 ) 

= 1 — 2zcos0„ + z 2 — 2(— l) n zzi sin0„, 
D 1 =D 1 {z,z 1 ) 

= 2izsin0„ - i(— l) n z 1 (l + 2zcos0„ + z 2 ). (25) 

This leads to the solution 

Xi 1 - X M 



Bm — B\ - + 

x¥ — x^ 1 



A+ - A_ 



(022-81 - Ol2-Dl)- 



A M-1 _ A M-1 



Xj 



D 



M 



X, 



A_ 



(onDi - a 2 \B 1 )- 



A 



A_ 
■A J 



A_ 



(26) 



where A± = z(l — 2: 2 )e ± *^™' are the eigenvalues of the 
2x2 matrix in ( |23| ) . After some algebraic manipulation, 
this yields the expression (^) quoted in the preceding 
section. 



IV. THE (2M - 1) x N MOBIUS STRIP 

We consider a (2M— 1) x N Afobius strip in this section. 

In order to make use of results of the preceding sec- 
tions, we start from the 2M x TV strip of Sec. 2, and let 
spins in the two center rows of the strip (the Mth and 
(M + l)th rows) having interactions Kq = K/2. The ex- 
ample of a 4 x 5 lattice with these interactions is shown 
in Fig. 3. Then, by taking K\ — 00 {z\ = 1) as described 
in Sec. 1, this lattice reduces to the desired (2M — 1) x N 
Mobius strip of a uniform interaction K. 
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FIG. 3. Labelings of a 4 x 5 Mobius strip which reduces to 
a 3 x 5 Mobius strip upon taking K\ — 00. 

Following this procedure, we have 
Zf M \ N {K) = 2( 2M - 1 ) Ar (coshX) 4 ( M - 1 > JV 

x cosh 2Ar (K/2)G(z,z ,zi) , (27) 

zi = l 

where zq = tanh(if/2), and G(z, Zq, z\) is the generating 
function of closed polygons on the 2M x N Mobius net 
with edge weights as described in the above. 

The generating function G(z,Zq,Zi) can be evaluated 
as in the previous sections. In place of (|]), ( p^ ) and (pi]), 
we now have 



PL4(z, Zq, zi) = \J det.4(z, z , z±) 



2/V 



Y[ yJdetA M {z,z ,zi;< 



(28) 



with 

A M {z,z ,zr,<t>n) = 

/ B(z) a ,i{z) 

00,-1(2) B(z) 00,1(2) 

00,-1(2) 

V 



B(z) 00,1(2) 
00,-1(2) C(z 0> zi)J 

(29) 
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Then fph becomes 



G(z,z ,zi) = - (l - i)PiA(z,z ,izi) 



+ (l + i)PiA(z,zo,-izi) 



(30) 



The evaluation of det Am (z, z , z\ \ <p n ) can again be done 
recursively. Define as before Bm — det Am (z, z , Z\\ <f> n ) 
and Dm the determinant of Am with the fourth row and 
column removed, one obtains again the recursion rela- 



tions ( p3[ ) and arrives at precisely the same solution (26) 
but now with a different initial condition 



Bi = Bi(z Q , zi), 



Di = L»i(z ,zi), 



(31) 



where the functions B\ and D\ are defined in d25|) . After 
some algebra, this leads to 



2yAM-l)N 

c\ sinh Mt(<p n ) — c 2 sinh(Af — l)t(4> n ) 



PfA(z,z ,Zi) = [z(l-Z 

N 

><n 

where 



sinht( 



(32) 



Cl 



1 - z 



cos0„ + (— l) n Zi sin</>. 
c 2 = 2z -> 1 + z cos0„ + (— l) n z\ sin 




(33) 



The substitution of (|3^) into ( |30| ) and (|27]) now completes 
the evaluation of the partition function for a (2M— 1) x N 
Mobius strip. 



V. THE KLEIN BOTTLE 

The Ising model on a Klein bottle can be considered 
similarly. We consider first a 2 M x N lattice £, con- 
structed by connecting the upper and lower edges of the 
Mobius strip of Fig. 1 in a periodic fashion with N extra 
vertical edges. As in the case of the Mobius strip, it is 
convenient to let the extra edges have interactions A 2 . 
The solution for a uniform interaction A is obtained at 
the end by setting K\ = K 2 = K. 

The Ising partition function for the Klein bottle now 
assumes the form 

Z«£ N (K,K U K 2 ) = 2 2MN (co S hK) 4MN - 2N 

x (cosh K-l cosh K 2 ) N G Kln (z , z x , z 2 ) , (34) 



where 



G Kln (z, Zl ,z 2 )= 



(35) 



generates all closed polygons on the 2M x N lattice £ 
with edge weights z — tanh K, z\ — tanh K\ , and z 2 = 
tanh A' 2 . The desired partition function is given by 



Z$g tN (K, A, A) = 2 ZMN (cosh A) 4Miv G Kln (z, z, z). 

(36) 

Again, it is necessary to first write G Kln (z, zi, z 2 ) as a 
multinomial in z, zi, z 2 in the form of 



G Kln (z,z!,z 2 )= J2 T m , n (z)z?z2, 



(37) 



where T m ^ n (z) are polynomials in z with strictly positive 
coefficients. 

The evaluation of G Kln (z, z\, z 2 ) parallels that of 
G(z, zi) for the Mobius strip. One first maps the lat- 
tice C into a dimer lattice Co by expanding each site 
into a city of 4 sites as shown in Fig. 2. Orient all A 
and Ai edges of Co as shown, and orient all A 2 edges in 
the same (downward) direction as the K\ edges. Then 
this defines an 8MN x 8MN antisymmetric matrix ob- 
tained by adding an extra term to A(z, zi) given by (|l4|), 
namely, 



A K[n {z, z u z 2 ) = A(z, zx) + b(z 2 ) ® G'j 
Here, 



b(z 2 ) 



H 



2N- 



(38) 
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(39) 



Then, in place of theorem (||) , we now have 

N 

' Kln (z, zi, z 2 ) = Y] emenT , m,n(^)^r^: 



PfA 1 ^ 



2 J 



(40) 



from which one obtains in a similar manner the result 

G Kln (z, Zl ,z 2 ) = 

- [pM Kln (z, izi, -iz 2 ) + PfA Kln ' 



'(z, -izi,iz 2 ) 



closed polygons 



-iPfA Kln (z, izi, iz 2 ) + iPf A Kln (z, -izi, -iz 2 ) . (41) 

To evaluate the Pfaffian (|40|), we note that the matrix 
( p8| ) can again be diagonalized in the {2A} subspace, 
yielding 
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PfA Kln (z,z u z 2 ) = ^detA K ^(z,z u z 2 ) 



2N 



z,z\,z 2 \ 



(42) 



where 

Am u (z, z x ,z 2 ; (f> n ) = A M (z, zi;<fi n ) 

+i{-l) n+1 b{z 2 )®G' M . (43) 

Now we expand detA^ ln in z 2 . Since setting z 2 = the 
determinant is precisely Bm and the term linear in z 2l 
the {4,4} element of the determinant, is by definition 
Dm, one obtains 

deL4^ n (z, z x , z 2 ; cj> n ) = B M + i(-l) n z 2 D M , M > 2, 

(44) 



where Bm and Dm have been computed in (|26j). This 
leads to 



PfA Kln (z,z 1 ,z 2 ) = (1 

N 

><n 

n=l 

where 

c(z, Zl, z 2 ) = 



Z\Z 2 



N 



z(l-z 2 ) 



MN 



sinh(M + l)t - c(z, Zi,z 2 ) sinhMi 



sinh t 



(45) 



1 



(l + z 2 )(z 4 + Zl z 2 ) 



z(l - z 2 ){z 2 + ZlZ 2 ) 

+2z(z 4 — z\z 2 ) cos (f> n 
+2(-l) n (z 1 + z 2 )z 3 sin 0„j. (46) 

Setting z\ = z 2 = z in ( pd| ) and using (ff5|), we obtain 
after some algebra 



G Kln (z,z,z) = z(l-z 2 ) 



MN 



N 



\ [ 2coshMt(0„) 



where 

D{4>n) 



( sinhM^n) " 

+Im n (isew^ 



(l-z 4 )-2z(l + z 2 )cos< 



(47) 



z(l-z 2 ) 



-4i(-l) n 2 sin^ n 



(48) 



and Im denotes the imaginary part. The substitution 
of ([£?]) into (H^) now completes the evaluation of the 
partition function for a 2M x iV Klein bottle. 
For the 2x2 Klein bottle, for example, we find 



PL4 Kln (z, Zl , z 2 ) = l + z 4 + 4( Zl + z 2 )z l - 2(z( + z 2 2 )z 2 
+2 Zl z 2 (l + z 2 ) 2 - 4ziz 2 (zi + z 2 )z 2 
+z 2 z 2 (l + z 4 ), 



G Kln (z, zi, z 2 ) = 1 + z 4 + 4(zi + z 2 )z 2 + 2(zf + z$)z 2 
+2 Zl z 2 (l + z 2 ) 2 + 4z 1 z 2 (z 1 + z 2 )z 2 
+z 2 1 z 2 (l + z i ), (49) 



which can be verified by explicit enumerations. 

For a (2 M -ljxfl Klein bottle we can proceed as 
before by first considering a 2M x N Klein bottle with 
interactions K,K\,K 2 and, within the center two rows, 
interactions K$ = K/2 as shown in Fig. 3. This is fol- 
lowed by taking K\ —t oo and K 2 — K. Thus, in place 
of (]34|), we have 

Zfrn-iA*) = 2 (2M " 1)A '(cosh^)( 4M - 3 ) JV 

x cosh 2A, (^/2)G Kln (^, z , 1, z), (50) 

where zq — tanh(if/2) and G(z, zo, Z\, z 2 ) generates 
polygonal configurations on the 2M x N lattice. Then, 
as in the above, we find 

G Kln (z,z ,zi,z 2 ) = 

1 [ P fA Kln (z, z , izi, -iz 2 ) + PL4 Kln (z, z , -iz x ,iz 2 ) 



-iPL4 Kln (z, zo, iz\,iz 2 ) + iPf A^ ln (z, z , -iz\, -iz 2 ) 

(51) 

where PL4 Kln (z, Z q, Z\, z 2 ) is found to be given by the 
right-hand side of J32]), but with 

ci = (1 + Zq)(1 - ziz 2 ) - 2z (l + ziz 2 ) cos<^)„ 
-2(-l)"(zi + z 2 )z sin</>„, 
1 



C2 



z(l-z 2 ) 



(z 2 + ziz 2 ) (1 - z z) 2 + (z - z )' 



+2(z - z )(l - z z) (z 2 - ziz 2 ) cos^) n 



+ (-l)"(z 2 zi + z 2 )sin< 



(52) 



expressions which are valid for arbitrary z, Zo, z x , and z 2 . 
For z = tanh(if/2), the case we are considering, (|52]) 
reduces to 



2z 



2z 



1 - ZlZ 2 - z(l + Z\Z 2 ) COS0„ 

-(-l) n z(z x + z 2 )sin0„ , 
z 2 + ziz 2 + z(z 2 - z\z 2 ) cos 0„ 
+ (-l)"z(z 2 zi + z 2 )sin</> Ir 



(53) 



which reduces further to (33) after setting z 2 = 0. The 
explicit expression for the partition function is now ob- 
tained by substituting (^Tj) into (|50|). 



VI. THE BULK LIMIT AND FINITE-SIZE 
CORRECTIONS 

In the thermodynamic limit, our solutions of the Ising 
partition function give rise to a bulk free energy 



G 



/bulk (AT) = lim 



1 



M.N- 



2MN 



In Z(K) 



(54) 



identical to that of the Onsager solution ||. Here, Z(K ) 
is any one of the 4 partition functions. Indeed, using the 
solution (111] ) for the 2M x TV Mobius strip, for example, 
one obtains 



A" 



/buik(^) = C{K) + lim (27V)- 1 V t(</> n ) 

n=l 

= C{K) + L f t (<j>) 

2tt 2 Jo Jo 
2cosh2i4:coth2i ; s:-2(cos6' + cos0) , (55) 



In 



where C(K) = [ln(2 sinh 2K)]/2. This expression is the 
Onsager solution (steps leading to the last line of ( p35| ) 
can be found in |ll|). It is well-known that /buik(AT) 
is singular at the critical point sinh2Ar c = 1 or 2K C = 
ln(V2 + 1). 

For large M and N, one can use the Euler-MacLaurin 
summation formula to evaluate corrections to the bulk 
free energy, an analysis first carried out by Ferdinand and 
Fisher (l2| for the Kaufman solution of the Ising model 
on a torus. Generally, for large M and N, we expect to 
have 

lnZ 2MiN (K) = 2MNf hulk (K) + Nc^K) 

+2Mc 2 ^,K) + c 3 (^,K) + --- (56) 

where £ = N/2M is the aspect ratio of the lattice. For 
the purpose of comparing with the conformal field theory 
[Q , it is of particular interest to analyze corrections at the 
critical point. Following |l^] as well as similar analyses 
for dimer systems p],|l3|, we have carried out such anal- 
yses for our solutions as well as for the solution of the 
Ising model under cylindrical boundary conditions HQ . 

For the 2M x N Mobius strip, for example, one starts 
with the explicit expression ( |Tl| ) of the partition function, 
and uses the Euler-MacLaurin formula to evaluate cor- 
rections to the bulk free energy. The analysis is lengthy, 
even at the critical point K c . We shall give details else- 
where and quote here only the results: 



_ Mob 



= 0, 



= 1 - K r = -0.087 618 



2<?1(oko 



■In 



1 + 



fl3(0|ig/2)-fl4(0|»g/2) 
2<? 3 (0|if) 



(57) 



'. sin (p + y 1 + sin 
= 0.353 068 ... 
and $i(u|T), i = 2,3, 4, are the Jacobi theta functions |i"4| ] 

oo 

M u \ T ) = 2 Y1 i (n ~^ )2 cos ( 2n - 

n=l 

oc 

i? 3 (u|r) = 1 + 2 ^ 9 ™ 2 cos 2nu, 

n=l 

oc 

i? 4 (tt|r) = l + 2^(-l) n 9 " 2 cos2nu, (58) 



n=l 



with q = e i7rT . For the 2M x N Klein bottle, we find 



c 2 (£,K c ) = 



In 



2^1(0|2^) 



i? 2 (0|2i^4(0|2^) 
1 



2t? 3 (0|2«0 



(59) 



If one takes the limit of N — ► oo (M — > oo) first while 
keeping M (A^) finite, one obtains 



lim — ln Z 2M ,n(K c ) = 2Mf bulk (K c ) + a 

N—>oc iv 

+A 1 /2M + 0(1/M 2 ), 

lim 7777 ln Z 2M,n{K c ) = Nf huik (K c ) + c 2 
M— >oo zikz 

+A 2 /iV + 0(l/7V 2 ), (60) 

where Ci, c 2 , Ai, A 2 are constants. We list our resluts in 
the table below. Also listed are values for the Ising model 
with toroidal boundary conditions taken from 12 , and 
values for the cylindrical boundary conditions computed 
using the solution of |?j . Our values of Ai and A 2 imply 
a central charge 



c = 1/2 



(61) 



for the Ising model. This is consistent with the conformal 
field theory prediction M. 



Mobius Klein Cylindrical Toroidal 



d c™ 5 " 



„Mob 





c 2 

Ai tt/48 tt/12 tt/48 tt/12 
A 2 tt/48 tt/48 tt/12 tt/12 



where 
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VII. SUMMARY 

We have solved and obtained close-form expressions 
for the partition function of an Ising model on finite 
Mobius strips and Klein bottles with a uniform interac- 
tion K. The solution assumes different forms depending 
on whether the width of the lattice is even or odd. For a 
2M x N Mobius strip, where 2M is its width, the parti- 
tion function Z^° N (K) is given by ( ]Tl| ) with F± given by 
(|l2|). For a (2M — 1) x N Mobius strip we employ a trick 
by first considering a 2M x N lattice and then "fusing" 
it into the desired lattice by coalescing two rows of spins. 
The resulting partition function Z^lfL 1 N {K) is given by 
( p7j ) in which the generating function G(z, Zq, z%) is ( |30| ) 
with the Pfaffians given by (32). 

For a 2M x N Klein bottle the partition function 
%2M-i n{K) i s given by (|3^ ) in which the generating 
function G Kln (z, z, z) is given by @. For a (2M-1) x N 
Klein bottle the partition function 2|a?-i n(K) i s given 
by ( |50| ) in which the generating function G Kln (z, Zq, 1, z), 
z a = tanh(if/2), is computed using (|5l]). All solutions 
yield the same Onsager bulk free energy (|55|). 

We have also presented results of finite-size analyses 
of our solutions as well as those of the Ising model un- 
der cylindrical and toroidal boundary conditions at crit- 
icality. The analyses yield a central charge c = 1/2 in 
agreement with the conformal field prediction Q]. 
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